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Abstract 

The searching for the stable patterns in the evolution of cellular automata is im- 
plemented using stochastic synchronization between the present structures of the 
system and its precedent configurations. For most of the known evolution rules with 
complex behavior a dynamic competition among all the possible stable patterns is 
established and no stationary regime is reached. For a particular rule coded by the 
decimal number 18, a self-synchronization phenomenon can be obtained, even when 
strong modifications to the synchronization method are applied. 
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1 Introduction 



Cellular automata (CA) are extended systems, discrete both in space and 
time. The simplest one dimensional version of a cellular automaton is formed 
by a lattice of N sites or cells, numbered by an index i = 1, . . . , A^. In each site 
a local variable cTj taking a binary value, either or 1, is asigned. The binary 
string a{t) formed by all sites values at time t represents a configuration of 
the system. The system evolves in time by the application of a rule A new 
configuration cr(t + 1) is obtained under the action of the rule $ on the state 
cr{t). Then, the evolution of the automata can be writen as 

a(t + l) = $[a(t)]. (1) 
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If coupling among nearest neighbors is used, the value of the site i, aiit + 1), 
at time t + 1 is a function of the value of the site itself at time t, o",j(t), and the 
values of its neighbors crj_i(t) and crj_|_i(t) at the same time. Then, the local 
evolution is expressed as 

ai{t + 1) = </>(cT,_i(i), ai{t), cTi+i(t)), (2) 

being a particular realization of the rule For such particular implemen- 
tation, there will be 2^ different local input configurations for each site and, 
for each one of them, a binary value can be assigned as output. Therefore 
there will be 2^ different rules 0. Each one of these rules produces a different 
dynamical evolution. In fact, the dynamical behaviors generated by the 256 
rules were already classified in four great groups. The reader interested in the 
details of such classification is addressed to the original reference where this 
development can be followed [1]. 

CA provide us with simple extended systems in which to prove different meth- 
ods of synchronization. A stochastic synchronization technique was introduced 
by Morelli and Zanette [2] that works in synchronizing two CA evolving under 
the same rule $. The two CA are started from different initial conditions and 
they are supposed to have partial knowledge about each other. In particular, 
the CA configurations, cr^(t) and cr^(t), are compared at each time step. Then, 
a fraction p of the total different sites are made equal (synchronized). The 
synchronization is stochastic since the location of the sites that are going to 
be equal is decided at random. Hence, the dynamics of the two coupled CA, 
a{t) = (cr^(t), cr^(t)), is driven by the successive application of two operators: 

(1) the deterministic operator given by the CA evolution rule $, — 
ma\t)ma\t)]\ B.nd 

(2) the stochastic operator Vp that produces the result rp[(T(t)], in such way 
that, if the sites are different {a} ^ af), then Tp sets both sites equal to 
al or erf with the same probability p/2. In any other case Tp leaves the 
sites unchanged. 

Therefore the temporal evolution of the system can be written as 

a{t + l) = {rpO^)[a{t)]^rpma\t)ma\tm. (3) 



A simple way to visualize the transition to synchrony can be done by displaying 
the evolution of the difference automaton (DA), 

S,{t)=\al{t)-am\- (4) 



2 



The mean density of active sites for the DA 



i=i 



(5) 



represents the Hamming distance between the automata and verifies < p < 
1. The automata will be synchronized when \im.t_^oo p{t) = 0. As it has been 
described in Ref. [2], two different dynamical regimes, controlled by the pa- 
rameter p, can be found in the system behavior: 

p < Pc—>^ lim^^oo p{t) 7^ (no synchronization), 
p > Pc^ lim(_»oo p{t) — (synchronization) , 

being p^ the parameter for which the transition to the synchrony occurs. When 
P ^ Pc complex structures can be observed in the DA time evolution. In Fig. 
1 typical cases of such behavior are shown near the synchronization transi- 
tion. Lateral panels represent both CA evolving in time and the central panel 
displays the evolution of the corresponding DA. When p comes close to the 
critical value p^ the evolution of 6{t) becomes more rarefied and reminds the 
problem of structures trying to percolate in the plane [3]. A method to de- 
tect this kind of transition, based in the calculation of a statistical measure of 
complexity for patterns, has been proposed in the Refs. [4]. 



2 Self-synchronization of cellular automata 

2.1 First self- synchronization method 

Let us now take a single cellular automaton [5]. If o"^(t) is the state of the 
automaton at time t, o-^{t) — o-{t), and cr^(i) is the state obtained from the 
application of the rule $ on that state, o"^(t) = $[(T(i)], then the operator Fp 
can be applied on the pair {cr^{t), cr^{t)), giving rise to the evolution law 



The update of the different sites between ^^{t) and o'^(t) takes place with 
probability p in a^{t). This updated array is the one that stays as a{t + 1). 

It is worth to observe that if the system is initialized with a stable config- 
uration for the rule $, ^[cr] = a, then this state a is also stable for the 
dynamic equation (6), and the evolution will produce a pattern constant in 
time. In particular, such behavior can be observed in rule 18 where the system 
evolves toward its only stable null pattern. However, in general, this stability 



a(t + 1) = T,[ia\t),a'm = r,[(a(t), <I.[a(t)])]. 



(6) 
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is marginal. A small modification of the initial condition gives rise to patterns 
variable in time. In fact, as the parameter p increases, a competition among 
the different marginally stable structures takes place. The dynamics makes the 
system to stay close to those states, although oscillating continuously and ran- 
domly among them. Hence, a complex spatio-temporal behavior is obtained. 
Some of these patterns can be seen in Fig. 2. 

2.2 A second self-synchronization method 

Now we introduce a new stochastic element in the application of the operator 
Tp. To differentiate from the previous case we call it Tp. The action of this 
operator consists in applying at each time the operator Fp, with p chosen at 
random in the interval (0,p). The evolution law of the automaton is in this 
case: 

ait + l)=t,[ia\t),a'm = r,Mt),mm]- (7) 
The DA density between the present state and the precedent one, defined as 
S{t) =1 a{t) — a{t — 1) I, is plotted as a function of p for different rules $ in 
Fig. 3. Only when the system becomes self-synchronized there will be a fall 
to zero in the DA density. Let us observe again that the behavior commented 
in the first self-synchronization method is repeated here. Rule 18 undergoes 
a phase transition for a critical value of p. For p greater than such critical 
value, the method is able to find the stable structure of the system. For the 
rest of the rules the freezing phase is not developed. The dynamics generates 
patterns where the different marginally stable structures randomly compete. 
This entails for the DA density to linearly decay with p (see Fig. 3). 

2.3 A third self-synchronization method 

At last we introduce another type of stochastic element in the application of 
the rule $. Given an integer number L, the surrounding of site i at each time 
step is redefined. A site ii is randomly chosen among the L neighbors of site i 
to the left, {i — L, . . . ,i — 1). Analogously, a site v is randomly chosen among 
the L neighbors of site i to the right, . . . , i+L). The rule $ is now applied 
on the site i using the triplet {ii,i,ir) instead of the usual nearest neighbors 
of the site. This new version of the rule is called ^l, being $l=i = $. Later 
the operator Fp acts in identical way as in the first method. Therefore, the 
dynamical evolution law is: 

a{t + l)^rp[{a\t),a'm^rp[{a{t),^4a{tm. (8) 

The DA density as a function of p is plotted in Fig. 4 for the rule 18 and in 
Fig. 5 for other rules. It can be observed again that the rule 18 is a singular 
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case that, even for different L, maintains the memory and continues to self- 
synchronize. It means that the influence of the rule is even more important 
than the randomness in the election of the surrounding sites. The system self- 
synchronizes and decays to the corresponding stable structure. Contrary, for 
the rest of the rules, the DA density decreases linearly with p even for L = 1 
as shown in Fig. 5. The systems oscillate randomly among their different 
marginally stable structures as in the previous methods. 



3 Conclusions 



Inspired in stochastic synchronization methods for CA, different schemes for 
self-synchronization of a single automaton has been proposed and analyzed in 
this work. Self-synchronization of a single automaton can be interpreted as a 
strategy for searching and controUing the stable structures of the system. In 
general, it has been found that a competition among all such structures is es- 
tablished, and the system ends up oscillating randomly among them. However, 
rule 18 is a singularity within this panorama, since, even with random election 
of the neighbors sites, the automaton is able to reach the stable configuration. 
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Rule 30 




Rule 90 Rule 110 




Fig. 1. Spatio-temporal patterns near the synchronization transition. Lateral panels 
represent both CA evolving in time and the central panel displays the evolution of 
the corresponding DA. Time goes from top to bottom. Lattice size is iV = 100 and 
the number of iterations is T = 250. The stochastic coupling between both CA is 
p = 0.23. 
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rule 90 




Fig. 2. Rule 90 has two stable patterns: one repeats the Oil string and the other one 
the 00 string. Such patterns are reached by the first self-synchronization method 
but there is a dynamical competition between them. In this case p = 0.9. Binary 
value is represented in white and 1 in black. Time goes from top to bottom. 
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Fig. 3. Mean density p vs. pmax = p for different rules evolving under the second 
synchronization method. The existence of a transition to a synchronized state can 
be clearly observed for rule 18. 
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Fig. 4. Mean density p vs. p for rule 18 evolving under the third self-synchronization 
method. The existence of a transition to a synchronized state can be observed despite 
of the randomness in the election of neighbors within a range L, up to L = 4. 
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Fig. 5. Mean density p vs. p for different rules evolving under the third 
self-synchronization method. The density of the system decreases with p. 
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